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Abstract 

We study the possible mixings between gauge vector fields and scalar 
fields through their self-energies, arising in models with two Higgs doublets. 
We derive the relevant set of Schwinger-Dyson equations and the Ward iden- 
tities that compel the longitudinal parts of the field propagators. Linear 
gauge is used and the results are given at all orders in perturbative theory, 
and some particular aspects of the one loop case are stressed. 
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1 Introduction 



The Standard Model (SM) is nowadays tremendously confirmed by many 
experimental results, even though there is not yet direct evidence for scalar 
Higgs particle. Alternative models like the two Higgs doublet model (THDM), 
especially the minimal supersymmetric standard model (MSSM) are exten- 
sively studied, experimentally and theoretically as well. It is now well settled 
that THDM suffer from large quantum effects, and one loop corrections (and 
beyond) are really important to improve the knowledge of such models (see 
for exemple 1 and references therein). Particle mixings through quantum 
loops are common features of THDM and SM and we will study some of 
them. 

The mixing between the photon and the Z fields is well known and is 
inherent to the Glashow -Salam - Weinberg (GSW) model. More precisely , 
there is not only a mixing between the transverse parts of the propagators 
of the photon and the Z but also a mixing between the longitudinal parts 
of the photon and the Z propagators and the G° neutral Goldstone boson 
propagator. Moreover a mixing between the charged W ± fields and the 
charged Goldstone G ± exists as well. 

Actually these mixings come from the SU[2]*U[1] structure of the gauge 
group and the Higgs mechanism, so in any model with this gauge group and 
a symmetry breaking "a la Higgs" one gets such mixings. In the GSW model, 
the mixings were studied thoroughly by Baulieu and Coquereaux 2 . The 
aim of this paper is to extend their study to the THDM, where the Higgs 
scalar sector is richer. Since we focus on the mixings coming from the gauge 
group and the Higgs mechanism , our derivation is also valid for the MSSM, 
without any change. Actually , we will examine two mixings. The neutral set 
contains the photon , the Z , the Goldstone G and the CP odd neutral Higgs 
scalar A . In the charged set , there are the gauge field W ± , its associated 
Goldstone G ± and the Higgs field H ± . 

In order to keep the results as general as possible , we will use the BRS 
symmetry to find the constraints on the two point Green's functions. Indeed 
the BRS transformations of the fields are due to the gauge structure and 
the gauge fixing term , which turn out to be the same in all the models we 
consider. For the same reason, we carry out a derivation valid at all orders 
in perturbation theory, that gives as a by-product the one loop results and 
shows the physical differences between first order and all orders. 

For the sake of generality, because there is a large variety of possible sce- 
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narii, we don't at all study the renormalization of the two point functions, 
mixed or not mixed. Indeed, one can choose different renormalization condi- 
tions, one can renormalize before or after breaking of the SU[2]*U[1] gauge 
symmetry, one can choose a renormalization scheme where the gauge vector 
mixing ( photon-Z ) is present or absent, and , in the case of the MSSM, 
the dimensional reduction could be more adapted than the dimensional reg- 
ularization. Moreover, if one wishes that the constraints we derive using the 
BRS symmetry be respected before and after the renormalization procedure, 
it will be necessary to choose it in such a way the gauge fixing lagrangian 
be renormalization invariant, that it is always possible 3 . Finally, we will 
only assume that the Green functions we use are regularized with a gauge 
invariant method. 

This paper contains three sections. The first one is mainly devoted to 
display the tools we will use in the other sections. We begin by our notations 
that we wish natural and compact in such a way that the Schwinger -Dyson 
equations are comprehensive and clear. We also recall some details on the 
BRS transformations, useful to derive the constraints induced by the BRS 
symmetry on the various two point Green functions. 

In the second section we study the charged case, where we have to take 
into account three fields, six two point functions and one constraint. It is in 
some sense a short introduction to the neutral case worked out in the third 
section, a much more complex case since we have to handle four fields , ten 
two point functions and three constraints. 

2 Preliminaries 

2.1 Quantum numbers and mixings. 

In the THDM (or MSSM), the C and P quantum numbers of the Higgs 
or Goldstone scalar fields and the vector fields are well known 4 . They are 
summarized in the following table: 
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The second column concerns the bosonic interactions when charge con- 
jugation and spacial parity are separately conserved, while the third column 
includes the fermionic sector which is CP invariant. In the case of charged 
fields, C is not a good quantum number, so it has been discarded. If a neutral 
vector field has the 1 PC quantum number assignments , then its deriva- 
tive d^V 1 is a neutral scalar with the 0~ PC assignments. That allows some 
mixings between scalars and vectors via their longitudinal parts but also it 
forbids the mixing between the CP even Higgs scalar fields and the photon 
or the Z. This is the true reason why the Standard Model Higgs scalar field 
H ( which has the same quantum numbers as the THDM ho and H ) does 
not couple to the neutral vector gauge fields. In any case, one finds that 
the neutral mixing set contains the four first particles of this table and the 
charged set involves all the charged particles. 



2.2 Notations 

We now define our notations for the two point functions. They are slightly 
different from the notations of reference 2 . 

First we denote by L and T the longitudinal and transverse projectors : 

= k»k"/k? andT"" = - IT , 

in such a way that the free propagator of a gauge vector field reads 

D i»> = -i(T» v /(k 2 - m 2 ) + f L"7(Jfc 2 - £m 2 )) 

where £ is the gauge fixing parameter. To be complete, the free propagator of 
a massive scalar field is D = i/(k 2 — m 2 ) while the propagator of a Goldstone 
field is D = i/(k 2 -£m 2 ). 
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We will note the complete two point Green function and the one 
particule ireducible (1PI) two point Green function, where i and j are field 
indices: 



G lJ =<0\ l (x)j(y)\0> . 

In our definition of the 1PI two point Green functions, the tadpoles are 
allowed since we cope with expressions before renormalization. We can now 
introduce the transverse and longitudinal parts of these Green functions ( 
in momentum space ), where M is a totally arbitrary mass dimensioned 
parameter: 

G% = -i{G%T^ + Gf^IT) ; Pt- = +i{Pp^ + I* IT) 



G% = - iWG L JM ■ T% = +ik»P l L j /M 



In the previous formulae , two Lorentz indices mean a vector-vector mix- 
ing, one Lorentz index means a vector-scalar mixing and no Lorentz index 
means a scalar-scalar mixing. With these formulae, the G functions ( like 
the D functions ) and P functions have respectively the canonical dimension 
-2 and + 2. 



2.3 The Schwinger - Dyson equations 

In pure QED the Schwinger -Dyson equation can be proved analytically, or 
diagramatically by "sticking" the product PD to the geometrical expansion 
of G in terms of D and P. In both cases one gets the well known result: 
G = D + DPG , whose diagrammatic form is: 

iwvwww^/vwwvw = mmm + m/vwvw(J]])m^^ 

This pure QED equation can be easily extended when mixings arise. Using 
a matrix formalism, one readily finds the general equation G = D + DPG, 
where G , D and P are now matrices. If a free field propagator is generically 
painted as a solid line, and a complete ( 1PI ) Green function is painted as 



■5 



a black ( grid ) disk between two solide lines , then the pictorial form of this 
matrix equation is as follows: 



+ 



As it is possible to get vector- vector , vector-scalar and scalar-scalar mixings 
, we can display three equations , where v (s) is the set of all the possible 
vector (scalar) fields: 



'i Z / Vi ViV m V m oi 



G Sk s t - DsJh + E D s k p s k v m G v m s l + J2 D s k Ps k s n G Sn s l ■ 



Introducing the transverse and longitudinal Green functions, one obtains 
four coupled equations, one tranverse and three longitudinal: 

G T ViVj = (Sij + E Pv,v m G T VmVj )/(k 2 - m 2 ) (2.1) 

men 



G L ViVj = (<% + E Pv.vAv, + k 2 E PisG L SnV jM")/(e - £m 2 ) (2.2) 
= £(E ^v m G^ m5l + E ^k^)/(A; 2 - £m 2 ) (2.3) 



met) nes 



G L SkSl = + k 2 E ^ G WM 2 + E ^^/(fc 2 - m\) (2.4) 

In the last equation, if Sk is a Goldstone boson, k 2 — ml must be changed in 
k 2 — im\. Of course, G and P are symmetric functions in their arguments. 
These four equations are the most general Schwinger -Dyson equations for 
mixings between spin zero and spin one particles. Their solutions give the 
complete two point Green functions G in terms of the 1PI Green functions 
P. 
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2.4 Ward identities. 



Since the THDM is a gauge model, Green functions are not independent and 
they are constrained by the BRS symmetry , the quantum level symmetry 
that encompasses the gauge symmetry valid at the classical level. Using 
the BRS transformations, one gets identities between the complete Green 
functions G, which in turn give identities between the 1PI Green functions 
P, when solved the above Schwinger -Dyson equations. 

First, one considers the following two point Green function : 

<0\ V *(x)F j (y)\0>, 

where Fj is the gauge fixing term associated to the gauge vector field Vj, 
and rji is the Faddev-Popov ghost associated to the gauge vector field Vj. 
By fermionic antisymmetry , this Green function is null, and so its BRS 
transform. Since the BRS transform of an anti-ghost rj* is the gauge fixing 
term Fj and the BRS transform of a gauge fixing term is proportional to the 
equation of motion of the corresponding anti-ghost , one gets 2,3 ' 5 

< OlF^Fjiy^O >= SijS(x-y) (2.5) 

This equation gives us the first set of Ward identities between the two 
point Green functions G . Since they only mix gauge vector fields and Gold- 
stone fields, these identities are valid in the GSW model and in any model 
with a larger Higgs sector as well. To be precise we will use in this paper the 
definitions of Bohm et al. 5 for the gauge vector fields and the gauge fixing 
terms, i.e.: 

= -sin9 w W^ + cosd w B^ = cos9 w W^ + sin6 w B^ ; 
W± = (WlTiWl)/V2 ■ 
F p = d^/y/S; F z = (d^ - im z G )/Ji ; 
F ± ^(d^W^iCm w G ± )/^. 
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A second set of Ward identities can be obtained with the BRS variation 
of the following Green function: < 0\r)*(x)Sj(y)\0 > when Sj is a Higgs scalar 
field. It is also a null function and its BRS variation reads: 

< 0)^(^)^(^)10 > - < 0\r ] *(x)S BR sS j (y)\0 >= (2.6) 

Such a formula cannot give simple relations between the G functions studied 
here, since the second term on the left side is a composite operator with 
an external ghost field. Because the fermionic ghost number is conserved, 
only this rj* can give a pole in the second term. Hence, multiplying by the 
inverse propagator of the scalar field Sj and putting it on its mass shell ( 
both operations noticed "S.o.s" ), the composite operator vanishes, and we 
just get identities between the G Green functions we consider here: 

<0\F, l (x)S J (y)\0> s .o,= (2.7) 

Moreover, when Fi is the photon or the Z gauge fixing term and Sj is the 
CP-odd neutral Higgs A , equation (1.4) gives us an useful result at one loop 
order 6 . Indeed the BRS transform of the A contains fields ( scalars and 
F.P ghosts ) in such a way that no tree level coupling exists with the r/ 7 and 
r] z ghosts . Hence at one loop order ( "o.l.o" index ) we get the following 
equations : 

[k 2 Gph]o.i.o =< 0\d^(x)A o (y)\0 > oXo = (2.8) 

and 

[k 2 G L zh + tiMm z G L gy ] . Lo =< 0\(d fl Z^x)-^m z G (x))A (y)\0 > . Lo = (2.9) 

which are true independently of any mass or energy prescription. 

It is to be noted that this second set of Ward identities (2.6), (2.7), (2.8) 
and (2.9) is of interest only for the extensions of the GSW model, as the 
THDM or MSSM, where the Higgs sector is rich enough to contains a neutral 
CP odd A and a charged scalar H ± , together with correct quantum numbers. 
It is also to be stressed that the identities due to the BRS symmetry are not in 
general renormalization invariant; but the way we derive the identities (2.5), 
(2.7), (2.8) and (2.9) assures they can be preserved by renormalization, as 
we said in the introduction. 
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3 Charged fields 



The charged gauge boson W ± can be mixed with the Goldstone G ± and the 
Higgs H ± , but only by its longitudinal part . For the transverse part , there 
is no mixing at all and the equation (2.1), which is very similar to the photon 
one in pure QED , reads , omitting the ± index : 



G 



T 

WW 



[l + PLG T ww )/(k 2 -m 



w) 1 



the solution is obvious: 



G lw = V(^ 2 



P T ). 

WW J 



Actually the interest is in the six coupled longitudinal equations coming 
from (2.2), (2.3) and (2.4) where we also have discarded the L index : 



G ww — £[1 + P ww G wg + k 2 (P wg G gw + P wh G hw )/M 2 ]/{k 2 — £m 



Gwg £\PwwG W g ~\~ PwgGgg ~\~ P w hG hg\ I (k £ 



/// 



w) 1 



G W h — i[P ww G w h + PwgGgh + P w hGhh\/{k 2 — im 2 w ) ; 



Ggg = "[I + k 2 Pg W G W g/M 2 + Pg h G h g + P gg G gg ]/(k^ ~ £ 



Gg h = -[k 2 Pg W G wh /M 2 + PggGg h + F ^G ^] / '(^ ~ ^ ) J 



G hh = -[1 + k 2 P hw G wh /M 2 + P hg G gh + P hh G hh ]/(k 2 



A solution of this linear system is rather easy to find with an analytic soft- 
ware but the solution suffers potentially large simplifications and it is highly 
preferable to "prepare" the system as follows : 
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with the notations : 

an = k 2 - £m 2 w - £P WW ; a 12 = k 2 a 23 /M 2 ; a 14 = k 2 a 25 /M 2 ; 

«23 = —£,Pwg ', <225 = —£,Pwh ', 

^33 = -(k 2 - £m 2 w + P gg ) ; a 35 = -f-P fl/l ; 

aee = ~£(k 2 -m 2 h + P hh ) . 

Now it is easy to obtain the most compact solution of the longitudinal charged 
mixing: 

G ww = iM 2 [-P 2 gh + {k 2 - iml + P gg )(k 2 -m 2 h + P hh )]/D, 
G wg = iM 2 \P gh P wh - P wg (k 2 -m 2 h + P hh )]/D, 
G gg = -[£k 2 P 2 h + M 2 (k 2 -m 2 h + P hh )(k 2 - iml - iP ww )]/D, 
G wh = £M 2 [P gh P wg - P wh {k 2 - im 2 w - £P WW )]/D, 
G gh = \ik 2 P wg P wh + M 2 P gh (k 2 - Zm 2 w - ZP WW )]/D, 

G hh = -[ik 2 P 2 wg + M 2 {k 2 - im 2 w + P gg ){k 2 - im 2 w - £P WW )]/D, (3.1) 
where the denominator D is : 

D = M 2 P 2 gh {k 2 -im 2 w -iP ww ) 

+tk 2 \P 2 wg {k 2 -m 2 h + P hh ) + P 2 h (k 2 - {ml + P gg ) - 2P gh P wg P wh ] 



+ M 2 (k 2 - iml + P gg ){k 2 -m\ + P hh ){k 2 - im 2 w - ^P ww ). (3.2) 
At one loop order, it remains : 

Gww = i{k 2 - im 2 w + iP ww )/{k 2 - imlf, 

G wg = —£,P wg /(k 2 — ^m^) 2 , 

G g9 = -(k 2 -Zml-P gg )/(k 2 -Zm 2 J 2 , 

G wh = -£P wh /[(k 2 - m 2 )(k 2 - {ml)], 

G 9h = P gh /[(k 2 -m 2 h )(k 2 -Zm 2 J], 

G hh = -{k 2 -m\- P hh )/{k 2 - ml) 2 . (3.3) 

As a consequence of this one loop result, we see that the first three Green 
functions G ww , G wg and G gg become independent of P g h, P w h and Phh ■ 
Therefore, within the GSW model, these three longitudinal Green functions 
(at one loop order) have exactly the same form as we have got. 
The relevant Ward identity (2.5) in the charged sector is: 

<0\F ± (x)F ± (y)\0>=5 ±± 5(x-y) 
and reads in the momentum space : 

k 2 G ww ± 2im w k 2 G wg jM + e™ 2 w G gg = f (3.4) 

the ± sign depending on the field W ± . When we replace in this equation 
the above solutions, we get the following constraint between the various P 
Green functions: 

2k 2 P gh P wg P wh - k 2 P gg P 2 h - (Mm w P gh t k 2 P wh f - M 2 P ww P 2 h 
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+ (k 2 -m\ + P hh ) [M 2 {k 2 + P gg ){m 2 w + P ww ) - k 2 {P wg ± Mm w )] = 0. (3.5) 

It is possible to simplify the denominator in the solutions of the previous 
system thanks to the Ward identity; then we get the following formula : 

D = -k 2 {MP ghT Zm w P wh ) 2 

+M(k 2 -m 2 h + P hh )(M(k 2 - im 2 w f + Mk 2 P gg T 2£k 2 m w P wg + £ 2 Mm 2 w P ww ). 

Of course the Ward identity (3.4) will be automatically fulfiled if one uses 
such a denominator. 

Again , giving up the scalar Higgs H ± as in the GSW model , the Ward 
identity (3.5) is simpler and reads: 

k 2 {Mm w ± P wg ) 2 = M 2 (k 2 + P gg )[m 2 w + P ww ) (3.6) 

This equation is quite similar to the equation 3.12 in 2 . It is worth noting 
this equation is almost what we get from the full equation at one loop order, 
since we just have to discard P 2 terms to end with: 

Mm 2 w P gg t 2k 2 m w P wg + k 2 MP gg = (3.7) 

So this constraint on the self-energies is at one loop order the same in both 
models, GSW or THDM. In other words , at one loop order, the charged 
Higgs H ± decouples in this constraint, which could be expected from the one 
loop solution of the system (formulae (3.3)). 

The second constraint we have in the charged sector follows directly from 
the equation (2.6) and we get the relation between longitudinal self -energies, 
with an on shell Higgs field: 

M{m\ - iml)[m\P wh ± Mm w P gh ) + 

MP gh {£Mm w P ww =f m hPw g ) + ITT'hPwhi^nT'wPwg T MP gg ) = 0. 
At one loop order the previous equations obviously reads 7 : 

m 2 h P wh ± Mm w P gh = 
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4 Neutral fields 



Likewise the third section , we begin by the transverse neutral case where 
the mixing between the photon and the neutral boson Z appears; we have 
three coupled equations ( 2.1 ), omitting the T label: 

Gpp (1 + PppGpp + PpzG pz) / k , 

Gpz {PppGpz + PpzG zz) / k , 

Gzz = (1 + P P zGpz + PzzG Z z)/(k 2 - m 2 z ). 
The solutions are ; 

G pp = (k 2 -ml- P zz )/[(k 2 - P pp )(k 2 -m 2 z - P zz ) - P 2 Z ]; 

Gpz = Pp Z /[(k 2 - P pp ){k 2 -ml- P zz ) - P p 2 J; 

Gzz = (k 2 - Pp P )/[{k 2 - P pp )(k 2 -ml- P zz ) - P 2 Z }. 

The really involved part is the longitudinal case where we have to deal 
with ten coupled equations ; indeed we have to take into account the photon 
, the neutral gauge vector boson Z, the Goldstone scalar Go and the CP 
odd scalar Higgs A , so we get 10 (= 4+3+2+1) possible G and P Green 
functions. To be short, we will not fully write these equations like in the 
third section , since they are easy to derive from the general Schwinger - 
Dyson relations ( 2.2 ). 

These complete longitudinal Green functions (omitting now the L label) 
verify the three following Ward identities derived from equation (2.5) : 

k 2 G P p = i (4.1) 
MG pz = iim z Gp g (4.2) 
k 2 MG zz - 2iim z k 2 G zg + i 2 m 2 z MG gg = £M (4.3) 
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and we will solve our problem using these constraints. If one try to solve this 
10x10 system by brute force one gets a huge result. Once again, we prefer to 
solve it in the following compact form : 
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where the coefficients an read as follows: 



a i4 = k 2 a 46 /M 2 ; a 17 = k 2 a 49 /M 2 ; a 35 = /c 2 a 56 /M 2 ; a 38 = k 2 a 59 /M 2 ; 
an = k 2 - £P PP ; a 12 = -£P P z] a 33 = k 2 - £m 2 - £P ZZ ; 

«46 = — £,Ppg] °49 = —£Pph', ^56 = —£,Pzg;', «59 = —£,Pzh', 

am = -£(k 2 - im\ + P gg ); a 69 = -£P flh ; a 00 = -£(& 2 - m 2 + P hh ). 



Since we have to deal with a 10x10 matrix, we naively expect that the 
determinant is a polynomial of degree 10 in a^, but if we smartly perform 
the derivation , we get solutions as: 

Gij = Nij/5, 

where Nij are polynomials of degree 3 while 5 is a polynomial of degree 4, 
which is much better . Now using the first constraint , i.e. the photon Ward 
identity ( 4.1 ) , N pp k 2 — £5, we can replace 5 by a polynomial of degree 3. 
Then we obtain solutions such that G pp = (,/k 2 and the other functions G^ 
= Tiij j k 2 d, where now d and Tiij are both polynomials of third degree . Of 
course we still have to rewrite everything in terms of the P^- functions to end 
up with the final result: 

Gpz = -eiM 2 P pz (k 2 -£m 2 z + P gg )(k 2 -m 2 +P h h)-k 2 P pg P zg (k 2 -m 2 +P hh ) + 
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k 2 P gh {P ph P zg + P pg P zh ) - k 2 P ph P zh {k 2 - im\ + P gg ) - M 2 P 2 gh P pz }/k 2 d. 

G zz = -m 2 P 2 pg {k 2 -m 2 h + P hh ) - ik 2 (2P gh P pg P ph - P 2 vh ik 2 - im\ + P gg ))- 

M 2 (k 2 - iP PP ){P 2 gh - {k 2 -m 2 h + P hh )(k 2 - im\ + P gg ))]/k 2 d. 

G P9 = mM 2 P pz {P zg {k 2 -m 2 h + P hh ) - P gh P zh ) - ik 2 {P ph P zg P zh - P pg P 2 zh )- 

M\k 2 - im\ - iP zz ){P gh P ph - P pg (k 2 -m 2 h + P hh ))]/k 2 d. 

G zg = -^M 2 P pz (P gh P ph - (k 2 -m 2 h + P hh )P pg - ik 2 P ph {P ph P zg - P pg P zh ) + 

M 2 (k 2 - iP pp ){P gh P zh - P zg {k 2 -m 2 h + P hh ))]/k 2 d. 

G gg = ieP pz (2k 2 P ph P zh - M 2 P pz {k 2 -m 2 h + P hh )) + 

& 2 {P 2 zh {k 2 - £P pp ) + P 2 ph {k 2 - im\ - £P ZZ )) + 

M 2 (k 2 - iP pp ){k 2 -m 2 h + P hh )(k 2 - im\ - iP zz )\lk 2 d. 

G ph = e[^ 29 (A; 2 (P P / l P 2S -P ra P 2 / l )-M 2 P 9A P p2 )+eM 2 P p2 P 2/l (A; 2 -em2+P 99 ) + 

M\k 2 - im\ - iP zz ){P ph {k 2 - im\ + P gg ) - P gh P pg )}/k 2 d. 

G zh = e[^ P9 (A; 2 (P P ,P 2 / l -P P / l P 2 J-M 2 P^P p2 )+eM 2 P p ,P p2 (A; 2 -em 2 +P ra )- 
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M 2 (k 2 - ZP pp )(P gh P zg - P zh {k 2 - im\ + P gg ))]/k 2 d. 

G gh = \e{M 2 P gh P 2 z - k 2 P pz (P ph P zg + P pg P zh )) - £k 2 (P zg P zh (k 2 - tP pp ) + 

P pg P ph (k 2 - im 2 z - iP zz )) - M 2 P gh (k 2 - iP pp ){k 2 - im 2 z - iP zz )]/k 2 d. 

G hh = ieP pz (2k 2 P pg P zg - M 2 P pz {k 2 - im\ + P gg )) + £k 2 (P 2 g (k 2 - iP pp ) + 

Ppg^k 2 — C m l ~~ Z,Pzz)) + M 2 {k 2 — iPpp) {k 2 — £m 2 + P gg ) (k 2 — ^m 2 — £P ZZ )]/ k 2 d. 
The denominator is : 

d = ik 2 \2P gh P zg P zh - P 2 g (k 2 -m 2 h + P hh ) - P 2 h (k 2 - im 2 z + P gg )} + 

M\k 2 - im\ - iP zz ){P 2 gh - (k 2 -m 2 h + P hh )(k 2 - im\ + P gg )). 

The two other constraints due to the two last Ward identities give two 
cubic polynomial equations. The first one ( 4.2 ) is quite simple and reads: 
Mn pz — i^m z n pg = 0. In terms of the , we get : 

046(059069 ~~ 000^56) + 049(056069 — 059066) + 

M 2 a 12 (ao a m - a 2 m )) + i£m z [k 2 a 59 (a 4e a 59 - a 49 a 56 ) + 

M 2 (a 00 (a l2 a 56 - a 33 a 46 ) + 059(033049 - ai 2 a 59 ))] = (4.4) 
while in terms of self-energies we obtain : 

-M 2 (k 2 - m 2 )(k 2 - £m 2 z )(MP pz + im z P pg )- 
M[MP pz ((MP gg + itm z P zg ){k 2 - ml) + MP hh {k 2 - im 2 z ))- 
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(k 2 - £m 2 z )P ph (k 2 P zh + iMm z P gh )- 



P pg (tMm z (P hh (k 2 - im 2 z ) - £P zz (k 2 - m{)) + k 2 P zg (k 2 - m 2 h ))] + 



M 2 P pz {P gh {MP gh + %im z P zh ) - P hh {MP gg + i£m x P za ))+ 



P ph (k 2 P zh (MP gg + %im z P zg ) - MP gh (k 2 P zg - i£Mm z P zz ))+ 



P pg {MP hh {k 2 P zg + %iMm z P zz ) - k 2 P zh {MP gh + iim z P zh )) = 0. (4.5) 

For the third Ward identity ( 4.3 ), the result is not so simple as ( 4.2 ); 
in terms of , we find 

066^49 — 2a4 6 a49a 6 9) + M 3 /c 2 an(ag 9 — a 00 a 6 Q) — 
Mk 2 (k 2 (a 00 al 6 + a m al 9 - 2a 56 a 59 a 69 ) + M 2 a^(a 2 m - a oa 66 )) + 
2i£m z k 2 [k 2 a A9 (a A9 a 56 - a 46 a 59 ) + 
M 2 (a 00 (a l2 a A6 - a u a 56 ) + a 69 (a u a 59 - a 12 a A9 ))] + 
£ 2 Mm 2 z k 2 (a 33 a 2 A9 - 2a 12 a A9 a 59 + 011059) + 

£ 2 M 3 m 2 z a 00 (a 2 12 - an 033) = 0. (4.6) 

We don't write it in terms of the P^ Green functions since we can improve 
it thanks to the two other Ward identities (4.1 and 4.2). We briefly present 
the method. We first solve the equation (4.4), i.e. we find eiy as functions of 
the other clm ; 059 and a^g are not suitable because they are irrationnal and 
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complex functions. Secondly , we substitute these in the equation ( 4.1) 
to get 7 equations , more or less complicated , and not independent. After 
simplifications by factors proportional to the determinant ( so different from 
zero ), only the simplest equations are interesting; they are due to aoo , ^33, 
<256 and clqq. It turns out that the first of them , written in functions of the 
various P functions , gives only quadratic and cubic terms in P , which is 
not really convenient for an eventual loop expansion. Among the three other 
possibilities , the best choice we can do is due to a^Q for two reasons: it gives 
the most compact form and at one loop order it is independent of the other 
results. Now we just have to combine the result given by the equation (4.6) 
with the null identity just derived to get in terms of the P functions: 

-M 2 (k 2 - m 2 h )(Mm 2 z P gg + 2ik 2 m z P zg + Mk 2 P zz ) 
M(M 2 m 2 z P 2 h + k 2 {k 2 - m 2 h )P 2 zg - iMm z P hh {2k 2 P zg - iMm z P gg ) + 
2%iMm z {k 2 - m 2 )(P pp P zg - P pg P pz + 2ik 2 Mm z P gh P zh + k 4 P 2 h - 
M 2 {{k 2 - m\)P gg + k 2 P hh )P zz ) - 2iim z {k 2 P ph {P ph P zg - P pg P zh ) + 

M (Phh(.PpgPpz PppPzg) ~\~ Pgh(PppPzh -fp/i-fpz))) - !"" 
M(k 2 (P z h(P gg P z h — PghPzg) + Pzg(PhhPzg ~ PghPzh)) + 

M 2 P zz (P 2 h - P gg P hh )) = 0. (4.7) 

The last relations between the P Green functions come from the equation 
( 2.7) ; from the relation concerning the photon , one easily gets: 

M 2 (m 2 h - im 2 z )[{m 2 h - £m 2 g )P ph + P gg P ph - P gh P pg + S(P px P zh - P ph Pzz)] 

+i[M 2 (P g h(PpgPzz — PpzPzg) + Pgg(Pp Z Pzh ~ PphPzz)) 
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+ m 2 h P zg {P ph P zg - P pg P zh )\ = (4.8) 
and from the relation with the Z boson, one obtains: 

M 2 m\(m\ - tml)(m\P zh - iMm z P gh ) + 
M[i(m 2 h - £m 2 z ){m 2 h P ph {MP pz - im z P pg ) - MP pp [m\P zh - iMm z P gh )) 
+m 2 h {m 2 h {MP gg - %im z P zg )P zh - MP gh {m 2 h P zg - i^Mm z P zz ))\ 
Pzh{Ppg{fn 2 hPpg i£Mm z P pz ) - MP pp (MP gg - i£m z P zg )) 
+M 2 P gh {P p MMm z P pz - m\P gg ) + P pp [m\P zg - i^Mm z P zz )) 

+ m 2 P ph (P pz (M 2 P gg - %iMm z P zg ) - P pg (m 2 P zg - i£Mm z P zz ))} = (4.9) 

both equations at the energy of the A Higgs mass. In the Landau gauge ( 
£ = with our notations ), the two previous constraints turn out to be very 
simple without cubic terms in P^ and just reduce in : 

Pph( m h + Pgg) = PghPpg, 

and 

P zh (m 2 h + P gg ) = Pg h (P zg - iMm z ). 

The first order expansion of all these results ( solutions and constraints ) 
are easy to derive. The solutions of the linear system become : 

Gp P = £/k , 
G pz = ePpJ{k 2 {k 2 -im 2 z )) , 
G zz = £(k 2 - im\ + iP zz )/{k 2 - im\) 2 , 
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G P9 = -iP P9 /{k 2 {k 2 - {ml)) , 



G zg = -iP zg /{k 2 -im 2 z f , 
G gg = (-k 2 + ^m 2 z + P gg )/(k 2 -m 2 z ) 2 , 

G ph = -ZP ph /(k 2 (k 2 -m 2 h )) , 
G zh = -iP zh /{{k 2 -m 2 h ){k 2 -im 2 z )) , 
G gh = P gh /{{k 2 -m 2 h ){k 2 -im 2 z )) , 

G hh = (-k 2 +m 2 h + P hh )/(k 2 - m 2 h ) 2 , 

while the constraints derived from equations 4.1, 4.2 and 4.3 are simplified 
as follows 2 : 

P =0 
MP pz + im z P pg = 

Mm\P gg + 2ik 2 m z P zg + k 2 MP zz = 0. 

At last, the constraints 2.8 and 2.9 read ( without any energy prescription 
as we have proved in section 1-4 ): 

Pph — 0, 

and 

k 2 P zh + im z MP gh = 0. 

The most striking difference between general and one loop results is P pp = 
0, which means the self -energy of the photon is only transverse at one loop 
order. The second point to be stressed is a general feature: everything is 
mixed beyond one loop, the CP odd Higgs included. 
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5 Conclusion 



In this study we have derived in a very general way the Schwinger - Dyson 
equations for mixings between scalar and vector fields in 577(2) * U{1) gauge 
models like the SM or the MSSM. We have found the solutions of these equa- 
tions and given some constraints induced by the BRS symmetry they have 
to verify . Altough this study is rather technical, it can be interesting for 
several reasons. For example, many authors use the perturbation theory at 
more than one loop and can find useful informations in our work. This is 
also true when one uses renormalization group methods to improve first order 
results. An another interest one can find in this paper is to check codes used 
to generate Feynman diagrams and to calculate amplitudes. Indeed, since we 
only have considered two point functions, all the formulae we have derived 
at one loop order can be written in terms of BO Veltman functions, favouring 
analytic verifications. 
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